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ABSTRACT 

We discuss how the solar occultations of bright sources of energetic gamma rays can 
be used to extract non-trivial physical and astrophysical information, including the 
angular size of the image when it is significantly smaller than the experiment's angular 
resolution. We analyze the EGRET data and discuss prospects for other instruments. 
The Fermi Gamma Ray Space Telescope will be able to constrain the size of a possible 
halo around 3C 279 from observations it makes on the 8th of October each year. 

Key words: gamma-rays: theory. 



1 INTRODUCTION 

The brightest source in the sky almost at any wavelength, 
the Sun is very weak in high-energy [E > 100 MeV) gamma 
rays. This property can be used to study solar occultations 
of gamma- ray sources. 

The width of the point-spread function (PSF) of tele- 
scopes detecting photons at these energies is quite large, 
of order several degrees. The enormous exposure of the 
Fermi Gamma Ray Space Telescope (the telescope previ- 
ously known as GLAST) would partially compensate for the 
poor resolution; however, it would be almost impossible to 
directly measure the angular size of the image which may be 
smaller than the PSF width. On the other hand, energetic 
gamma-ray images of distant sources may indeed have a sig- 
nificant angular size due to the cascading of photons on the 
background radiation and magnetic deflections of the cas- 
cade electrons and positrons. It has long been known that 
one can obtain the angular size of stars from lunar occulta- 
tions, we suggest that one may determine the image size of 
gamma-ray sources screened by the Sun0. 

The current collection of known energetic gamma-ray 
point sources is scarce (~ 300 sources detected by EGRET), 
so only a few are expected to be on the strip on the sky such 
that they are screened by the Sun. Fortunately, the bright- 
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^ Interestingly, we note that the Moon is much brighter than 
the Sun in this energy band because of seconda ry emission from 
cosmi c rays hitting t he lunar surface, see e.g. iThompson et al] 
lll997l' ): lBrigidal ||2009^ . 



est EGRET source identified with an extragalactic object, 
3C 279, has an ecliptic latitude of 0.2° and is screened by 
the Sun on the 8th of October each year. It is 3C 279 which 
is the main subject of our discussion because, as we will see 
in Sec.[2l it represents a perfect target for this kind of study. 

The simplest and most direct effect of an extended im- 
age size would be the detection of fiux from the source dur- 
ing occultation. Such a result could also be the signal of 
the transparency of the Sun t o gamma rays possible in sev- 
eral scenarios of new physics (|Fairbairn. Rashba fc Troitskvl 
l2[)07i" ) : however the parameter space of particular models is 
strongly disfavoured by results of other experiments (e.g. 
lAndriamonie et al.l (|2007l ')'). 

In Sec. |3] we review our analysis of archival EGRET 
data of the 1991 occultation of 3C 279 during which a non- 
zero flux was indeed observed, although at a very low sta- 
tistical significance. 

With a sensitive enough telescope, a more detailed 
study of the light curve during ingress and egress would 
be possible. In Sec. U we discuss the potential of Fermi for 
this kind of a study and mention sources other than 3C 279 
while Sec. O summarizes our conclusions. 



2 POSSIBLE MECHANISMS FOR EXTENDED 
EMISSION 

Very high energy photons interact with other photons in the 
source, with photons along the path between the source and 
the earth and with photons near to the sun. These inter- 
actions result in the production of electrons and positrons 
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which themselves consequently interact with the gamma-ray 
background leading to the developement of electromagnetic 
cascades which result in a gradual decrease of the average 
energies of propagating photons. Ambient magnetic fields 
deflect the electrons and positrons and, as a result, the image 
of the source seen in gamma rays becomes extended. This 
kind of extended image may be observed for distant sources 
which emit very energetic photons (both the Universe and 
typical source environments are transparent for gamma rays 
below ~ 1 GeV) . Recently, th e MAGIC collaboration has re- 
ported (|Teshima et al.ll2007^ the detection oi E > 200 GeV 
emission from 3C 279. Given the expected absorbtion on 
the cosmic infrared background, this corresponds to an ex- 
trem ely high luminosity of the quasar at very high ener- 
gies l| Albert et al.ll2008l). M oreover, according to both the 
Hillas criterion l|Hillaslll984 ) and to bounds on the source 
param eters from energy losses (see e.g. IPtitsvna fc Troitskvl 
l|2008h l. the jets of 3C 279 might provide the necessary con- 
ditions for the acceleration of cosmic rays (CR) up to ultra- 
high energies (UHE), E > 10^^ eV. 3C 279 therefore seems 
to be an ideal candidate to search for an extended halo due 
to the production of secondary photons. 

Let us summarize briefly a few different scenarios which 
could lead to the formation of an extended image: - 

Inverse Compton effect in the source envi- 
ronment (< GeV). It has been po inted out some 
time ago llh aronian. Coppi fc Vo"i5ll994 ) that a halo of 
(sub-)GeV inverse-Compton photons may form around the 
so urce of very high energy gamm a rays. The estimates 
of lAharonian. Coppi fc V^ (|l994 ) for the case of 3C 279 
give rise to an expected angular size of the halo of ~ 0.2° 
and a flux from the halo approximately equal to the flux 
coming from the central point source. 

Synchrotron halo of UHECR sources (~GeV). 
For reasonable values of the magnetic field (> 10~^ G) 
around a source of UHECR, synchrotron photons con- 
tribute to a halo of angular si ze of a fraction of degree 
potentially detectable b y Fermi (|Gabici fc Aharonian|[2007l : 
lAtovan fc Dermeij|2008l ). 

Electromagnetic cascades on the intergalactic 
magnetic fields (~ 10'^ GeV). Energetic photons from 
distant sources undergo electromagnetic cascades when scat- 
tering off extragalactic backgro und light, resulting in ex - 
tended images for TeV sources (|Neronov fc Sem ikoz"2007^. 
Even stronger cascading is exp ected for secondary photons 
from UHECR sources (see e.g. iFerrigno. Blasi fcPe Marcol 
l|2005l) ): because of their higher energies, they can scatter 
on CMB photons as well as the infrared background, and 
the number density of CMB photons is much higher. How- 
ever, extended emission of this kind is too weak at GeV 
energies to lead to a detectable effect. 

Electromagnetic cascades on the solar radiation 
(~ 10^ GeV). The solar radiation is strongly concentrated 
in the optical band, corresponding to the thermal emission 
of 5800 K, that is usun ~ 1 eV. The pair production thresh- 
old is determined hy E = mg/cjsun ~ 260 GeV. The optical 
depth of the solar radiation with respect to the pair pro- 
duction for 260 GeV photons tangent to the solar surface is 
r ~ 0.1 and any secondary electrons and positrons produced 
in this way would be isotropised by the ~ 1 G magnetic field 
of the sun. Because of this, the extended halo as viewed from 
earth would be too weak to be observable (one may hope to 
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Figure 1. Normalized light curves of an extended source oc- 
culted by the Sun, assuming Gaussian model for the extended 
image. Different values of the radial source variance a are taken 
in units of degree. The time corresponds to the angle (in degree) 
between source center and the point of its minimal separation 
from the center of the solar disk, for the minimal separation of 
0.20°. The shadow region represents the occultation of the point- 
like source. 

detect such a halo when the Sun passes in front of the regions 
with significant diffuse emission at 260 GeV if such regions 
exist, but without the timing signature, this emission would 
be hard to detect). 

Formation of the halo is a random process and therefore 
one expects that the variability timescale of the halo should 
be determined by its physical size (light minutes for the 
solar neighbourhood but millions of light years for all other 
scenarios). 

We see that one might expect a halo around 3C 279 
in the energy band detectable by EGRET and Fermi (above 
100 MeV). This halo is of angular size ~ 0.1° while the point- 
spread function (PSF) of these instruments extends for sev- 
eral degrees. However, the halo could be detected when the 
Sun screens the bright central source. Study of the shape 
and spectrum of this halo would help to distinguish between 
various scenarios of its formation and therefore contribute 
to our understanding of the source engine, the source en- 
vironment and the extragalactic background radiation. The 
extended halo would reveal itself in a smooth falling of the 
flux when the Sun approaches the source and in a non-zero 
flux from the source while it is screened by the Sun. In Fig.[Tl 
we present examples of the lightcurves for various source ex- 
tensions. 



3 EGRET OBSERVATIONS 

In 1991, the solar occultation of 3C 279 occured within the 
fleld of view of EGRET. We have analysed publicly available 
EGRET data to test the conjecture that the flux from the 
quasar w as non-zero when th e source was screened by the 
Sun (Fai rbairn. Rashba fc Troitskv 20071 ); here we present 
more details and discussions related to the study. 

Given the coordinates of 3C 279 taken from 
the NASA/IPAC Extragalactic Database (NED) 
( http://nedwww.ipac . caltech. edu| ) and the coor- 
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dinates of the Sun calculated with the program 
PLANEPH IChapront fc Francoul 1 19981 '). we are able 
to establish that the source was screened by the Sun for 
8 hours and 34 minutes. The minimal separation between 
the quasar and the center of the solar disk was 0.20° while 
the angular radius of the Sun at that per iod during the 
earth's orbit was 0.2675° (|Kulikovskvll2004l 'l. At that time 
(viewing period 11.0), EGRET was pointed in the direction 
of 3C 273, some 5° away from 3C 279. The quasar was in 
a moderate state and was firmly detected in ga mma rays 
during that viewing period l|Hartman et al.lll999l ). 

The actual distribution of observed photons 
with their coordinates, energies and arrival times, 
as well as a record of the instrument's viewing 
modes, are available from the EGRET ftp site 
ff tp: //cossc.gsfc. nasa.gov/compton7data/egret ). 
To calculate the exposure map for the peri od of occultation , 
we made use of the EGRET software (jBlackburnI 1 19951 : 
iLaubenthal et aLlliggj ). 

To determine the source flux, one compares the dis- 
tribution of arrival directions of detected gamma-rays with 
the sum of background and point-source fluxes using a par- 
ticular model for the former and the instrument point- 
spread function (PSF) for the latter, both convolved with 
time- and direction-dependent experimental exposure. We 
denote the diffuse background flux of gamma rays from 
a given direction (a. 5) as B(a, 5). This flux was deter- 
mined bv lCillis fc HartmanI (|2005h from the analysis of the 
EGRET data. From B{a,5), we calculate -B(V'), the ex- 
pected flux of background photons within the angular dis- 
tance tj) from the source we study (it is the integral over 
the corresponding circle on the sky, of B{a,5) weighted by 
exposure ) . The PSF of EGRET, obtained by its careful cal- 
ibration l|Thompson et af] 1 19931 ). is energy dependent; one 
has to assume some spectral index a for the source. We use 
the PSF for E > 100 MeV and a = 2 ( the EGRET-measure d 
spect ral index of 3C279 is 1.96 ±0.04 (jHartman et al.lll999l ')') 
from ICiUis fc HartmanI (|2005h and denote it as p{ip). 

The number n^tj}) of observed photons from the circle 
of angular radius t/j centered at the source (a ladder-like 
function of a single variable '(p) is fitted by the sum of a 
slowly varying background and a sharp PSF, 

n(V') = hSi-iP) + NspW + s4^), 

i 

where b and Ns are parameters of the fit (determined by the 
standard le ast-squares m e thod. Following the original 3EG 
procedure, iMattox et al.l ()l996[ ). we keep the coefficient 6 
free motivated by possible temporal variations of the back- 
ground) and Si are contributions of known nearby sources 
(for which we use the 3EG values for this viewing period). 
The best-flt number of source photons Ns divided by the 
exposure determines the source flux. In the same way as the 
authors of the 3EG catalog, we select events with energies 
E > 100 MeV, TASC in coincidence at 6 MeV and distance 
to the source not exceeding 15° . These cuts correspond also 
to the PSF we use. 

Our approach is very similar to that used for the con- 
struction of the 3EG catalog (see iMattox et all (|l996l ) for 
more details of the method) apart from some details: 

(i) We use the updated maps of EG RET- observed diffuse 
background (|Cillis fc HartmanI |2005| ) while iHartman et al.l 
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Figure 2. The normalized light curve (blue line) for the Gaus- 
sian extended source with the extension providing the best fit 
to the EGRET data on 3C 279 occultation in 1991. The time is 
determined as in Fig. [T] The horizontal band represents the 68% 
confidence-level EGRET ofi'-occultation flux determined from the 
two-week period. Black (thin) error bars give the average flux 
during the occultation seen by EGRET; red (thick) error bars 
are expected for a Fermi observation, in the pointing mode, of 
one occultation with the same parameters. 



|1999|) used a theoretical model for the distribution of the 
background flux. 

(ii) The catalog construction used the counts distribution 
binned in two celestial coordinates; we use the unbinned 
distribution in one coordinate - dis tance from th e sour ce 
(both approaches were discussed bv IMattox et al.l l|l996l )). 
This is more appropriate in the case of small number of 
observed events. 

(iii) The 3EG procedure substracts all 3-sigma sources 
(in the catalog, they list only 4-sigma ones). We substract 
only the sources listed in the catalog, even if they are fainter 
than 4-sigma for this particular viewing period (in this way, 
we potentially exclude some sources which did not pass the 
4-sigma cut for this viewing period, or for any other viewing 
period or their sum, but were brighter than 3-sigma in this 
period). 

We performed the flt for the occultation period and 
for the rest of the viewing period. The best fits give the 
signal of Ns — 4.82 photons for the occultation and A^'^ — 
284.5 for the rest of the viewing period. Within the 68% 
containment angle of the PSF, the numbers of signal and 
background photons are roughly equal. Therefore we found 
some weak evidence for a non-zero point-source flux from 
the location of 3C 279 during the occultation, (6.212'!/) • 
10~^ cm~'^ s~^ , to be compared with the value obtained from 
the analysis of the rest of the same viewing period, (8.6 ± 
0.5) ■ 10"'^ cm~^ s"'^ (the value quoted in the 3EG catalog 
for this period is (7.94 ± 0.75) • 10"^ cm"^ s"^). The zero 
flux (point-like source) is excluded at « 98% confldence level 
(CL); no upper limit can be placed on the source extension 
because unsuppressed flux is well within the 68% CL error 
bars. Assuming a Gaussian extended image, the best fit for 
the radial extension is cr ~ 0.283°. The error bars in Fig. [2] 
give an idea of the size of statistical uncertainties both of 
this result and of potential future studies. 
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3EG 
name 



other 
name 



(TS) 



1/2 



1219- 


-1520 




13.2° 


<0.89 


0.0 


12294 


-0210 


3C 273 


10.5° 


<0.95 


0.4 


1230- 


-0247 


1229-021 


7.0° 


1.13±0.43 


2.9 


1234- 


-1318 




9.3° 


<0.89 


0.0 


12354 


-0233 




9.9° 


1.24±0.41 


3.5 


12364 


-0457 




11.9° 


<0.90 


0.3 


1246- 


-0651 


1243-072 


2.5° 


1.29±0.54 


2.7 


1310- 


-0517 




3.6° 


1.05±0.51 


2.2 


1339- 


-1419 


1334-127 


13.6° 


<1.08 


0.0 


1255- 


-0549 


3C 279 


0.0° 


7.94±0.75 


15.1 



Table 1. Potential confusing sources: 3EG sources witliin 15° 
of 3C 279. e is the angular offset from 30 279; F-r is the flux 
during v. p. 11.0 in 10~^[cm~^ ( TS)^^'^ is the significa nce of 

detection in the v. p. 11.0. Data from lHartman et al ] l|l999l ). 



Here we mention some possible subtleties which one 
should be aware of when trying to understand this result 
or similar future observations. 

Other nearby 3EG sources. The numbers quoted 
above were obtained with the subtraction of the expected 
contribution of a single nearby source detected at 3 sigma 
in the v.p. 11.0 (3EG J1235+0233, see the list of nearby 
sources in Table O. The result can in principle be con- 
fused by contribution of other sources. As a test, we changed 
the 3-sigma threshold adopted by EGRET to both 2 sigma 
and to 4 sigma without any significant change in the result. 
The procedure described above assumed that the flux of the 
sources being subtracted was constant during the viewing 
period. Clearly, an extreme flare of 3EG J1246— 0651, or 
3EG J1310— 0517, or both, exactly at the occultation time, 
could explain our result without 3C 279. Note however that 
on the time scale between one viewing period and another, 
these sources do not demonstrate significant variability so 
such a flare seems unlikely. 

The Sun. The solar surface could be a gamma-ray 
source due to its interaction with cosmic rays. Early EGRET 
studies put a 95% C.L. upper limit of 2.0 ■ 10~^ cm~^ s~^ on 
the flux of the quiet Sun (|Thompson et al" ]|l997h . A marginal 
detection of solar disk flux ~ 4 ■ 10~ ^ cm~^ s~^ has been 
reported bv lOrlando fc Strond (|2008l ). The theoretical ex- 
pectation of the fl ux of the disk of the quiet Sun is about 
2 ■ 10"* cm-2 s"^ (|Orlando. Retry fc Stron3l2007D . 

The extended em is sion of the Sun 
IIMoskalenko. Porter fc Digell l2006l : lOrlando fc Strond 
l2007l . l2008t ) cannot explain the observed excess: within 
the 68% containment width of the EGRET PSF the 
expected flux from the extended solar emission is about 
assuming the model for the solar 



obtained by lOrlando fc Strond (|2008l ) in the 10-degree 
circle around the Sun. We would like to note that the recent 
unpublished Fermi observations have been used to extract 
the solar flux from the 10° regio n around the Sun which 
is (4.59 ± 0.89) ■ 10"^ cm"^ s~^ (lBrigidall2009l ). in a good 
egreement with lOrlando fc Strond i 2008l 'l. This flux is low 
compared to that of 3C 379 and may be neglected within 
our current poor precision. 



Start, UT 


End, UT 


EGRET photons 






(separation from the Sun) 


20'U4™ 


20*^15™ 


27.2°, 32.3°, 21.7°, 19.4°, 17.8° 


20'' 48™ 


20''49™ 




2l''59™ 


22''00™ 


10.8° 


22'' 07™ 


22''08™ 





Table 2. Solar flares during the occultation. 



of 



Solar flares. Solar flares are 
gamma rays; the BATSE 



sources 
records 



(http : //umbra. nascom.nasa.gov/batse/batse_years .html ) 
list four weak flares in the occultation time (see Table ^ . 
The photons detected by EGRET during these flares were 
separated from the Sun by at least 10.8°, so they most 
probably do not contribute to the point-source flux (the 
68% width of the EGRET PSF at E > 100 MeV is 3.3°). 

The Moon. One more gamma-ray source was nearby 
during the time of these EGRET observations: the Moon 
was in 6° to 9° from the Sun during the period of interest. 
The gamma-ray {E > 100 MeV) flux of t he Moon in 1991 
was (3 .6±0.9)-10~^ cm~^ s"^ (see Fig. 2 of lThompson et al.1 
([1993)). Given the separation, the PSF width and the flux 
we conclude that the lunar contribution cannot explain the 
observed excess. 

Possibility of misidentification. Though 3C 279 is 
considered as one of the best EGRET identifications, one 
still cannot exclude the possibility that the gamma-ray ex- 
cess is due to a source mis-identified with 3C 279. This 
would-be actual source, if located just 20 arc min away, 
would not be screened by the Sun. The best-fit position of 
the EGRET source associated with 3C 279 is indeed dis- 
placed from the position of the quasar as seen in other wave- 
lengths, but this position is evenly more deeply screened by 
the Sun during an occultation. 



4 PROSPECTS OF FUTURE OBSERVATIONS 

Let us estimate the ability of new experiments to observe 
the solar occultations of gamma-ray sources, notably that 
of 3C 279. AGILE cannot be pointed to the Sun because of 
configuration of its solar panels (M. Tavani, private com- 
munication). Atmospheric Cerenkov telescopes cannot be 
pointed at the Sun as they would be destroyed. The sen- 
sitivity of MILAGRO is insufficient to detect the source in 
8.5 hours. The Large Area Telescope (LAT) of Fermi may 
however be used to observe the occultation succesfuUy. 

In the survey mode, Fermi will scan the sky ro- 
tating continuously, so the short-period exposure to 
a given point in the sky is not too large. The sensi- 
tivity may be estimated using the Fermi web service 



(http:/ /www-glast. slac.stanford.edu/software/IS/glast_lat_performance. htm 
For the pointing mode, the average on-axis effective area 
for E > 100 MeV, weighted with the ~ E~^ spectrum, is 
~ 4700 cm^. Further precision may be gained by repeating 
the observation each year. 

To estimate the ability of the instrument to detect non- 
zero flux during the occultation, we have to assume a partic- 
ular value of the total flux of the source. 3C 279 is strongly 
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Figure 3. The expected number of source photons from 3C 279 
during the occultation, seen by Fermi, versus the radial extent of 
the Gaussian image, for 8.5— hour observations in the survey mode 
(full hne), pointing mode (dashed line), combination of three ob- 
servations in the pointing mode (dash— dot— dotted line) and for 
1— hour bin of the latter combination (dash-dotted line). The ver- 
tical line indicates the best-fit value of the extension from the 
EGRET data. The left and right scales are described in the text. 



name 


be 


3EG 


GEV 


VHE 


AX J1809. 8-2333 


-0.099 


J1809-2328 


yes 


no 


3C 279 


+0.186 


J1255-0549 


yes 


MAGIC 


W28=M20 (?) 


-0.020 


J1800-2338 


yes 


HESS 



Table 3. (EGRET and FERMI gamma-ray sources potentially 
eclipsed by the Sun. fee gives t he ecliptic latitude calculated from 
the FERMI bright source list jAbdo et al.ll2009l'l. "GEV" mean s 
E > 1 GeV detection by EGRET jLamb fc MacombI Il997l) . 
"VHE" means E > 100 GeV detection. 



variable; however the extended image should not be vari- 
able unless it is formed near the Earth and not near the 
quasar. While formation of the halo near the quasar is bet- 
ter motivated physically (cf. Sec.[2J, the opposite is in bet- 
ter agreement with the EGRET result (Sec. [Sjl because the 
non-variable flux of 3C 279, roughly estimated as a minimal 
flux over EGRET vie wing periods, is ~ 8 x 10~* cm~^s~^ 
l|Hartman et al.l 1 1999! ) , much lower than the best-fit flux 
during the occultation. Figure [3] shows the expected num- 
ber of photons from the source during the occultation as 
a function of the source extension for various Fermi ob- 
servations for both scenarios. The left-handed scale (AV) 
corresponds to the total flux of the Gaussian image of 
83.7 X 10~® cm~^ s~\ the average flux of 3C 2 79 over nine 
viewing periods l|Casandiian fc Grenieij 120081 '): the right- 
handed scale (NV) corresponds to the non-variable flux. The 
number of background photons within the 68% containment 
width of PSF is approximately equal to the number of source 
photons in the latter scenario. 

Two more EGRET sources with ecliptic latitudes be < 
0.25° (see Table [3]) had been classified as unidentified in the 
3EG catalog. Further studies suggested potenti al identifica- 
tions ; they also have been detected by FERMI (|Abdo et al.l 
I2OO9I ) The study of the solar occultation may help to de- 
termine their coordinates with higher precision, testing the 
identification. We expect that these (and maybe other) 



sources will become potential targets for angular-size mea- 
surements. 



5 CONCLUSION 

It will be interesting to try and measure the angular sizes of 
images of energetic gamma-ray sources by means of obser- 
vation of their solar occultations. The best target is 3C 279, 
whose occultation happens each year on the 8th of Octo- 
ber. EGRET observations made during such a period did 
not exclude the unsuppressed flux of the quasar when it was 
screened by the Sun. The sensitivity of the Fermi telescope 
is high enough that if the flux was unsuppressed during oc- 
cultation, it could be observed more definitively than with 
EGRET. Fermi can also constrain the angular size of the 
image even in the survey mode, and is capable of obtaining 
a light curve by the combination of several observations in 
the pointing mode. This would help to constrain models of 
particle acceleration and magnetic fields in and around the 
quasar. If the flux during the occultation exceeds the non- 
variable flux of the source (as it is slightly favoured by the 
EGRET data), it would mean that either the extended im- 
age is formed relatively nearby or the Sun is partially trans- 
parent for the point-like gamma-ray emission (both options 
would mean a discovery of some unconventional physical or 
astrophysical phenomenon). The same method may be ap- 
plied to refine the coordinates and/or to estimate the angu- 
lar size of images of other gamma-ray sources screened by 
the Sun. 
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